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Abstract 

Space-borne interferometric gravitational wave detectors, sensitive in the low-frequency (milli- 
hertz) band, will fly in the next decade. In these detectors the spacecraft-to-spacecraft light-travel- 

O , 

times will necessarily be unequal, time-varying, and (due to aberration) have different time delays 
on up- and down- links. Reduction of data from moving interferometric laser arrays in solar orbit 



will in fact encounter non-symmetric up- and downlink light time differences that are about 100 
times larger than has previously been recognized. The time-delay interferometry (TDI) technique 
uses knowledge of these delays to cancel the otherwise dominant laser phase noise and yields a vari- 
ety of data combinations sensitive to gravitational waves. Under the assumption that the (different) 



o 

■ up- and downlink time delays are constant, we derive the TDI expressions for those combinations 

CO ! that rely only on four inter-spacecraft phase measurements. We then turn to the general problem 

o 

"q ' that encompasses time-dependence of the light-travel times along the laser links. By introducing a 

cr 

set of non-commuting time-delay operators, we show that there exists a quite general procedure for 
deriving generalized TDI combinations that account for the effects of time-dependence of the arms. 
By applying our approach we are able to re-deriye the "hex-free" expression for the unequal-arm 
Michelson combinations X\, first presented in |1|], an d obtain the generalized expressions for the 
TDI combinations called Relay, Beacon, Monitor, and Symmetric Sagnac. 
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I. INTRODUCTION 

Future space-borne gravitational wave (GW) observatories, such as LISA [2], will have 
sensitivity in the low-frequency band and will use time-delay interferometry (TDI) to can- 
cel laser phase noise. All the original papers on TDI considered a configuration of three 
spacecraft interchanging coherent laser beams, and tacitly or explicitly assumed the array 
to be at rest in an inertial system. TDI was treated in Euclidean 3-space with a universal 
time, in which the velocity of light is c and isotropic. Recipes were given for combining 
data (time series) separately recorded at the various spacecraft, delayed by transit times 
calculated from the inter-spacecraft separations Lj (i = 1,2,3V, in order to remove the oth- 
erwise overwhelming phase noise of the laser sources . The aim is possible 
detection of incident gravitational waves of galactic or cosmic origin. 

The LISA mission [2] will have three spacecraft orbiting the Sun in a triangular array 
with the Li ~ 5 x 10 6 km, and GW detection capability in the band 10 -4 — 1 Hz. Several 
TDI Michelson-like and Sagnac-like reduced laser-phase-noise-free data streams will have 
different responses to secondary phase noise sources and to two polarizations of incoming 
gravitational waves from different directions. A recent study of a linear array of three 
spacecraft in a single solar orbit (SyZyGy) 10| uses a TDI combination sensitive to a single 
polarization of incident gravitational waves, and two others sensitive solely to secondary 
system noises. 



In an important development, Shaddock JJ| noticed that rotational motion of an array 



results in a difference of the light travel times in the two directions around a Sagnac circuit. 
Two time delays along each arm must be used, say L, and L i for clockwise or counterclockwise 
propagation as they enter in any of the TDI combinations. Shaddock emphasized the need for 
careful distinguishing of primed and unprimed delays in the TDI combinations for Michelson- 
like combinations, and, to eliminate laser noise from the Sagnac-type combinations when 
the array is moving, he presented new TDI variables related to those originally given by 
being "double differenced". 

Cornish and Hellings also considered the effect of rotation of the LISA triangle around 
its centroid on the TDI combinations, and reported the new data combinations. Summers 



12j further pointed out that the LISA array is not rigid, 



13j and Cornish and Hellings 
that Li and L i not only differ from one another but can be time dependent (they "flex"), 
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and that again the laser phase noise (at least with present laser stability requirements) can 
enter at a level above the secondary noises. For LISA, and assuming ~ lOm/sec jl^ . 
they estimated the magnitude of the remaining frequency fluctuations from the laser to be 
about 30 times larger than the level set by the secondary noise sources in the center of the 
frequency band. This may not be as serious a problem with SyZyGy jlo| . 

Finally Shaddock et al. jl| addressed the "flexing" complication by showing that it 
becomes of higher order if the sequence of various time delays in the new doubly differenced 
Sagnac combinations is respected in the TDI recipe, and they introduced a new doubly- 
differenced Michelson-type TDI combination to achieve the same result. They stressed that 
although all these combinations are considerably more complicated than those originally 
given for a non-moving array, and their GW response functions are similarly complex, the 
final sensitivity — calculated from GW signal strengths and secondary phase noises — is 
unaffected. 

All the analyses above, however, assumed the clocks onboard the three LISA spacecraft 
to be synchronized to each other in a reference frame attached to the LISA array. It is well 
known Q], however, that the spacetime geometry - here the Sagnac effect - prevents the self- 
consistent synchronization of a network of clocks by transmission of electromagnetic signals 
in a rotating reference frame. This implies that the time adopted by the LISA onboard 
clocks and used for TDI has to be referenced to an inertial reference frame and that the 
onboard LISA receivers have to properly convert time information received from Earth to 
the time in this inertial reference frame. Within this frame, which we can assume to be 
Solar System Barycentric (SSB), the differences between back-forth delay times that occur 
are in fact thousands of kilometers, very much larger than has been previously recognized 
by us or others. The problem is not rotation per se, but rather aberration due to motion 
and changes of orientation in the SSB frame. 

In Section II, we further discuss the need for synchronizing the LISA clocks with respect to 
a common inertial reference frame (SSB), and the resulting GW response transfer functions. 
We turn in Section III to the derivation of the four-link TDI combinations valid for constant 
time delays. We first obtain the "unequal-arm Michelson" response, X, as an example of 
how time-delay operators can be used for deriving TDI data combinations. The operator 
formalism for TDI was introduced by Dhurandhar et al. We use it in conjunction 

with the usual subscripted delay notation to achieve a systematic understanding of the 
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"Relay" (U,V,W), "Beacon" (P,Q,R), and "Monitor" (E, F, G) combinations. With laser 
stabilization at a level somewhat improved from that used in the original LISA study j^, 
these combinations, now involving different up- and down- link delays, will satisfy sensitivity 
requirements. 

In Section IV, however, we go on to use delay time/operator notation to derive "second 
generation" TDI combinations, which account for both the inequality and time dependence of 
the back/forth optical paths. Following Shaddock et al. the resulting doubly-differenced 
combinations, immune to first order shearing (flexing, or constant rate of change of delay 
times) are denoted Xi,Ui, Pi, Ei , i = 1,2,3. All these new combinations suppress the 
nominal LISA laser phase noise to levels lower than those of the secondary (proof-mass 
and optical-path) noise sources, and their gravitational wave sensitivities are the same as 
previously computed for the stationary case. For completeness, we calculate the remaining 
shearing effect on the doubly-differenced versions of the system-noise-monitoring combina- 
tion (, denoted (i, (2, (3. Laser noise enters these combinations multiplied by sm 2 (nfL), 
where / is a Fourier frequency in the LISA band 10~ 4 — 1 Hz. We plot the laser noise in Q 
for the nominal LISA system and show it as a result also to be below the level of secondary 
noises. 



II. ABERRATION, TIME DELAYS AND SYNCHRONIZATION OF THE LISA 
CLOCKS 

The kinematics of the LISA and SyZyGy orbits brings in the effects of motion at several 
orders of magnitude larger than any previous papers on TDI have addressed. The instanta- 
neous rotation axis of LISA, and the SyZyGy array, both swing about the Sun at 30 km/sec, 
and on any leg the transit times of light signals in opposing directions, say Lj and L\ (c = 1), 
can differ by as much as 1000 km. Aberration due to LISA's orbit about the Sun dominates 
its instantaneous rotation. This observation reinforces the requirement that the new TDI 
combinations of Section III and IV must be used. Indeed, Lj and L • interchange periodically 
and so are also time dependent; this effect is however of order 0.1 m/sec and is dominated 
by the effect of shearing ( "flexing" ) already recognized. 

This large motional effect has been overlooked because intuitively up/down laser links 
between two spacecraft moving inertially on parallel geodesies certainly appears symmetric 
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in a co-moving frame. The spacecraft are then seen "at rest" and the elapsed light times, 
or delays, in either direction are the same. Consider however an inertial frame in which two 
spacecraft are moving with speed V along a line, with constant separation. The times of 
transit of a photon from one to the other, forward or back, clearly must differ by 2VL/c. 
This is just an extreme case of aberration. There is no paradox! We have taken the speed 
of light to be c and isotropic in both frames, and special relativity has taught us that that is 
fine so long as we properly re-synchronize the clocks that we use as time coordinates when 
changing frame (using light beams!) The spatial and temporal separations of two successive 
events along the null world line of a ray of light depend on choice of frame. Rays traveling 
in opposite senses between two moving spacecraft yield different separations in all frames 
except the co-moving one. 

An orbiting array is best described not by attempting a sequence of co-moving tangent 
"rest" frames, but rather in the barycentric non-rotating Euclidean frame moving with the 
Sun (of course we ignore tiny general relativistic distortions). The usual time coordinate 
of positional solar system astronomy in principle uses clocks such that c is isotropic. In 
LISA and SyZyGy data at the three spacecraft will undoubtedly be taken and time-tagged 
in this Solar System Barycentric frame, and all the up/down delay times used in the new 
TDI combinations must be calculated from the coordinates of emission and reception events 
in the SSB inertial frame. (This is exactly parallel to the time synchronization problem, 
and its resolution, that has been met by the designers of the GPS satellite array in 
geocentric orbit). 

Since the motion of the LISA array around the Sun introduces a difference between (and 
a time dependence in) the co-rotating and counter-rotating light travel times, the correct 
expressions for the GW contributions to the various first-generation TDI combinations will 
differ from the expressions valid for a stationary array (^j. The magnitude of the corrections 
introduced by the inequality of the light-travel times is proportional to the product between 
the time derivative of the GW amplitude and the differences between the actual light travel 
times. At one mHz, for instance, the correction to the expression of the signal valid for a 
stationary array is five orders of magnitude smaller. Since the amplitude of this correction 
scales linearly with the Fourier frequency, we can completely disregard this effect (and also 
the weaker effect due to the time dependence of the light travel times) over the entire LISA 
band. 
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FIG. 1: Schematic diagram of LISA configurations involving six laser beams. Optical path delays 
taken in the counter-clockwise sense are denoted with a prime, while unprimed delays are in the 
clockwise sense. See text for details. 



It is clear, however, that over many months of continuous observation of a quasi-periodic 
signal, the TDI responses have to account for the motion of the array around the Sun (and 
relative to the GW source), which introduces secular modulations in the phase, frequency, 
and amplitude of the GW responses jl7|. 



III. THE FOUR-LINK TDI COMBINATIONS: CONSTANT TIME DELAYS 

The notation we will adopt is the same as used in the paper by Shaddock et al. Qfl, (i.e 
it is different from the original TDI notation, e.g. Ref. |5j.) We distinguish time-of-flight 
delays by denoting with a prime those taken in the counter-clockwise sense and unprimed 
delays in the clockwise sense (see Fig. 1). 

There are six beams exchanged between the LISA spacecraft, together with the six phase 
measurements Sjj = 1,2,3) recorded when each transmitted beam is mixed with the 
laser light of the receiving optical bench. The phase fluctuations from the six lasers, which 
need to be canceled, can be represented by six random processes py, where is the phase 
of the laser in spacecraft j on the optical bench facing spacecraft i. In what follows we 
assume the center frequencies of the lasers are all equal, and denote it with Vq. Explicitly: 
S23 is the one-way phase shift measured at spacecraft 3, coming from spacecraft 2, along 
arm 1. The laser phase noise in s 2 3 is P32(t — Lf) — P2s{t), where we take c = 1, so that 
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L\ is the light time in the direction from spacecraft 2 to spacecraft 3. Similarly, S32 is the 
phase shift measured on arrival at spacecraft 2 along arm 1' of a signal transmitted from 
spacecraft 3. The laser phase noise in s 32 is p 23 (t — L\) — p 32 {t)i where L\ is the light time in 
the sense from 3 to 2 along arm 1 . For the further delays used in the TDI combinations we 
use the same conventions, being careful to distinguish light travel along arms with primes 
or not, depending on the sense of the measurement. For example, our notation for delaying 
the time series s 32 (t) by the clockwise light time in arm 1 would be 332,1 while delaying by 
the counterclockwise light time in arm l' would be S32.1'. As before, we denote six further 
data streams, Xy = 1,2,3), as the intra-spacecraft metrology data used to monitor 
the motion of the two optical benches and the relative phase fluctuations of the two lasers 
on each of the three spacecraft. The phase fluctuations of the lasers and optical benches 
enter into the measurements sy and in the following form (henceforth disregarding 
contributions from other noise sources and the gravitational wave signal) 



S31 
S21 

721 



Pl3 - f0 "2 ■ A 13 

P12 + von 3 - A 12 
P21 -psi-2 v h 3 ■ A 2 i + ii\ 

P31 - P21 + 2 U fl 2 ■ A31 + 



Psx + v h 2 - A 3 i 
P21 - ^0 n 3 ■ A 2 i 



(1) 
(2) 
(3) 
(4) 



In the above equations we have denoted with the phase fluctuations introduced by the 
optical fibers used for exchanging the laser beams between adjacent benches, and with the 
vector random processes Ay the phase fluctuations introduced by the mechanical vibrations 
of the optical benches. 

In order to simplify the derivation of the new TDI combinations, we note that by sub- 
tracting equation Q from Q we can rewrite the resulting expression (and those obtained 
from it by permutation of the spacecraft indices) in the following form 

1 

2 
1 

2 
1 



2 [T21 



T31 



- F32 - T12J 



T13 - T 23 



P31 + ^0 n 2 ■ A31 
P12 + ^0 n 3 ■ A 12 
P23 + vofi 1 - A 23 



P21 - vq n 3 • A 2 i 
P32 -vofti- A 32 

Pl3 -VQU2- A 13 



(5) 
(6) 
(7) 



If we now define the following combinations of laser and optical bench noises appearing in 
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equations (JTJI7J) |l| 

01 = Psi + tV A 3 i , (8) 
0i = P2i - ^3 • A 21 , (9) 

together with those obtained by permuting the spacecraft indices, it is possible to reduce 
the derivation of the new TDI combinations to the equivalent problem of removing the three 
random processes, 0i, <f) 2 , and 03, from the following six linear combinations of the one-way 
measurements sy and r^: 

Vn = s 2 i - - [r 32 - ri 2 ], 3 ' = 02,3' - 0i , V3i = s 3i + ^ [t 2 i - r 3 i] = 3j2 - 0i , (10) 

Vl2 = S12 + - [r 32 - T 12 ] = 0l i3 - 02 , ^32 = S 3 2 - ^ b~13 ~ T23], 1' = 03,1' ~ 02 , (H) 
»7l3 = S 13 - - [T 2 l - r 31 ] )2 ' = 01,2' - 03 , ^723 = ^23 + ^ [^13 ~ T 23 ] = 02,1 ~ 03 • (12) 

A. The Unequal- Arm Michelson 

Here we derive the unequal-arm Michelson combination, X, valid for the rigid-rotation 
case. We use X as an example for deriving TDI data combinations by using an alternative 
and powerful method based on the use of properly defined time-delay operators. 

The X combination relies on the four measurements r/i 2 , i] 2 i, ^13, and r/ 3 i. Note that the 
two combinations 7721 + 7/12,3', ^31 + 7/13,2, which represent the two synthesized two-way data 
measured onboard spacecraft 1, can be written in the following form 

7/21+7/12,3' = (D 3 ,D 3 -I) 01, (13) 
V31 + Vl3,2 = (D 2 D 2 > - I) 01 , (14) 

where we have denoted with Dj the time-delay operator that shifts by Lj the function it is 
applied to, and with I the identity operator. Note that in the stationary case any pairs of 
these operators commute, i.e. DiDj> — DyDi = (while they do not when the delays are 
functions of time l^J, jl]). 

From equations (|13 | I14 p it is easy to derive the following expression for X, by requiring the 
elimination of <f>i 

X = [D 2 D 2 , - I] fax + ?/i2,3') - [(Ds'Ds - I)} (r/31 + 7713,2) 

= [(7/31 + 7/13,2) + (t/21 + 7/i2,3'),2'2] - [(t/21 + 7/12,3') + (VSl + V 13, 2) ,33'] (15) 
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After replacing equations (fTUl ITTT IT2"j) into equation (fTBT) . we obtain the final expression for 
X valid in the case of rigid rotation of the LISA array [11[ 




As pointed out in [13j and equation (fT3j) shows that X is the difference of two sums of 
phase measurements, each corresponding to a specific light path from a laser onboard space- 
craft 1 having phase noise 0i. The first square-bracket term in equation (fTo^) represents a 
synthesized light-beam transmitted from spacecraft 1 and made to bounce once at space- 
craft 3 and 2 respectively. The second square-bracket term instead correspond to another 
beam also originating from the same laser, experiencing the same overall delay as the first 
beam, but bouncing off spacecraft 2 first and then spacecraft 3. When they are recombined 
they will cancel the laser phase fluctuations exactly, having both experienced the same total 
delays (assuming stationary spacecraft). 

B. The Relay 

The TDI "Relay" configurations were called (U, V, W) (equation (A4) of P). In what 
follows, let us consider, as a specific example, the U combination, which has to rely only 
on the four measurements 7731, 7712, 7732 and 7723. The idea we will follow for identifying the 
expression for U is to select combinations of some of these four measurements that contain 
only one phase noise. By then applying iteratively the time-delay procedure we introduced 
for the X combination, we will be able to remove all the phase noises fa,i = 1,2,3. Note 
that the obvious combinations that contain only one of the three phase noises fa are the 
synthesized two-way Doppler data measured onboard spacecraft 2 and 3. They in fact 
contain only the phase noises 02 and 03 respectively. Since the remaining two measurements 
7712 and 7731 can be combined in such a way as to eliminate the phase noise fa, we can start 
with the following set of three data combinations 



7?32,1 + ??23 = [DlDv ~ I) 03 , 
V23,V + V32 = [DyDx - I] fa ■ 



Vl2 + 7731,3 



D 3 D 2 fa - 02 , 



(17) 
(18) 
(19) 
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It is then easy to see that the expression for U is given by the following linear combination 
of the properly delayed equations (fTTl fTSl [TUJ) 

U = [DyDx - I] (t] 12 + 7731,3) + (7723,1' + ^32) - D 3 D 2 (7732,1 + V23) , 

= (7/12,11' + 7/31,311') - (t/12 + 7/31,3) + (7/23,1' + 7732) - (7732,123 + 7723,23) , (20) 

which, in terms of the one-way measurements sy and 7y, becomes 

^ — s 31,311' — $31,3 + Si2,ll' — S12 + S23,l' + s 32 — 532,123 — 523,23 

+ -^[(t 2 i - T3l),3U' - (T21 - r 3l),3 ~ (l~32 ~ T12) + (?32 - T 12 ) ,1V 

+ (Tl3 - T 2 3),1'123 - (Tl3 — T 23 ) ^ (21) 

with V, W obtained by cycling the spacecraft indices. 
C. The Beacon 

In the "Beacon" combination, one spacecraft transmits (only) to the other two while those 
other two exchange one-way beams as usual. These were called the (P,QM combtaations, 
depending on which spacecraft was the transmit-only element of the array |5j]. In order to 
derive the expression for P, which involves only the four data streams 1712, 7713, r/32, and 7723, 
we will proceed according to the above considerations, and use in this case the following 
data combinations 

7/12,2' - 7/13,3 = £>303 - Ai'02 , (22) 
7/32,1 + 7/23 = [DxDy - I] 3 , (23) 
7/23,1' + 7/32 = [D V D 1 - I] 02 • (24) 

By taking advantage of the commutativity of the delay operators in this constant time delay 
case, it is easy to see that the expression for P is given by the following linear combination 
of the properly delayed equations ((221 EO 121) 

P = D 3 (7732,1 + 7723) - D 2 > (7/23,1' + 7/32) - [Di'Dx - I] (7712,2' - 7713,3) , 

= (7/32,13 + 7723,3) - (7/23,l'2' + 7/32,2') ~ (t/12,2'11' - 7713,311') + (t/12,2' ~ 7/l 3 ,3) • (25) 
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Equation (|25|) can be rewritten in terms of the one-way measurements s^, r 2 



P 



Sl2,2' — s 13,3 — ^32,2' + S23,3 + 532,13 ~ S 23,1'2' + 513,311' — s 12,2'll' 

+ -[(^"21 - 7"3l),2'3 - (T21 - T31) ,11'2'3 + (r 32 - Ti 2 ),2' ~ (^32 - Tl2),ll'2' 

+ (Tl3 - T23),3 - (Tl3 ~ ^Xll's] > 



(26) 



with Q, R obtained by cycling the spacecraft indices in Eq. (|2*6^l . 
D. The Monitor 

Similarly, there are three combinations where one spacecraft is listen-only Q]. In order 
to derive these "Monitor" combinations (E,F,G) (equation (Al) of [a]), let us consider the 
following combinations of the four data streams that enter into E 



Similarly to the derivations made for the two previous combinations, it is easy to see that 
the expression for E is given by the following linear combination of the properly delayed 
equations (E3 EH 123) 



7723,1' + 7732 = [DyDi - I]4> 2 ■ 



?732,1 + 7723 = [DtDy - I) 03 , 



7731 - 7721 



D 2 (f)3 ~ As' 02 



(27) 
(28) 
(29) 



E = D 2 (7732,1 + 7723) - D 3 > (7723,1' + 7732) - [DyDi - I] (7731 - 7721) , 

= (7732,12 + 7723,2) - (7723,1'3' + 7732,3') - (7/31,11' - 7/21, l'l) + (t/31 - 7/2l) , 



(30) 



which, in terms of the one-way measurements Sy and becomes 



E 



$32,12 + S 23,2 — S23,1'3' — £32,3/ + S21,l'l — s 31,ll' ~~ s 21 + $31 

+tt[(t"21 - 7-31) - (t 2 i - Tzi),iv + (r 32 - ri 2 ),3' - (r 32 - ri 2 ),ii'3' 




(31) 



with F, G obtained by cycling the indices. 



E. The ( Combinations 



In all the above, we have used the same symbol (e.g., X for the unequal-arm Michelson 
combination) for both the rotating (i.e. constant delay times) and stationary cases. This 
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emphasizes that, for these TDI combinations, the forms of the equations do not change going 
from systems at rest to the moving or rotating case. One need only distinguish between the 
time-of-flight variations in the clockwise and counter-clockwise senses (primed and unprimed 
delays) . 

In the case of an array at rest there is one symmetric data combination that cancels 
exactly all laser noise and optical bench motions and has the property that each of the 77^ 
enters exactly once and is lagged by exactly one of the one-way light times. We called this ( 
equation (3.5)) and showed how to take advantage of its relative immunity to GWs in 
order to assess on-orbit instrumental noise performance and distinguish instrumental noise 
from a confusion-limited background Although now the rotation of the array breaks the 
symmetry and therefore the uniqueness of a "£-like" combination, it has been shown (jll|. 



12|) that there still exist three generalized TDI laser-noise- free data combinations that have 
properties very similar to (, and which can be used for the same scientific purposes. Here 
we derive these combinations, which we call (Ci, C2, (3), by applying our time-delay operator 
approach. As we will see in the following section, our derivation will automatically identify 
the "correct" order of the delays that has to be applied to the one-way data. In other words, 
the expressions lead to an order of time delays such that even with shearing the remaining 
laser noise is below the level identified by the secondary noise sources. (1 will not have to 
be further generalized. 

Let us consider the following combination of the r]ij measurements 

7/13,3' - 7/23,3' + 7721,1 = [D&Dy - £>l] 01 , (32) 
7731,1' - 7/32,2 + 7/12,2 = [D^D 2 - D V ] 01 , (33) 

where we have used the commutativity property of the delay operators in order to cancel 
the </>2 and 03 terms. Since both sides of the two equations above contain only the 0i noise, 
(i is found by the following expression 

Ci = {D 3 ,D 2 > - £>i] (7731,1/ - r/32,2 + 7712,2) - [D 2 D 3 - D v ] (7713,3/ - 7723,3/ + 7721,1) . (34) 



12 



In terms of the one-way measurements and Ty, equation becomes 

Cl = [ s 31,l' — s 32,2 + Si2,2],2'3' _ [ s 13,3' — s 23,3' + s 21,l],32 
— [ s 31,l' — s 32,2 + s 12,2],l + [ s 13,3' ~~ s 23,3' + s 21,l],l' 

+ 2(( T 32 — Tl2),22'3' — ( r 32 ~ 7~12),21 + ( r 32 ~ 7l2),13'32 _ (t~32 ~ 7"l2),13'l' 

+ {Tl3 ~ T23),22'3'l' ~ (^13 — 7"23),211' + ( r 13 ~ r 23),233' — ( r 13 — T23),3'l' 

+ (T21 - r 3 i) ]2 2'33' - (721 — T3l),ll'] (35) 

to geth er with its cyclic permutations. (This expression for (i was given (but not derived) 
in and independently by Q|.) If the light-times in the arms are equal in the clockwise 
and counterclockwise senses (e.g. no rotation) there is no distinction between primed and 
unprimed delay times. In this case, £i is related to our original symmetric Sagnac ( by 
Ci — C23 — C,i- Thus for the practical LISA case (arm length difference < 2%), the SNR of 
(i will be the same as the SNR of (. 



IV. THE SECOND-GENERATION TDI COMBINATIONS 

Generalizations of the original unequal-arm Michelson, (X,Y,Z), and Sagnac, (a,/?, 7) 
TDI combinations to an array with systematic spacecraft velocities, showing that they ef- 
fectively cancel all laser phase noises, have been derived in Here we complete that set 
of TDI combinations by deriving generalized expressions for the "Relay", "Beacon", and 
"Monitor" combinations that are unaffected by the rotation and time-dependence of the 
light-path delays. These TDI combinations rely only on four of the six possible one-way 
measurements LISA will make, and for this reason they add robustness and trade-off op- 
tions to the LISA design. Like the unequal-arm Michelson combination X\ 1], these new 
combinations involve the four one-way inter-spacecraft measurements at 16 different times. 

The order of the time-delay operators now becomes important for laser phase terms. 
The operators can no longer be permuted freely to show cancellation of laser noises in 
the TDI combinations (they no longer commute!). In order to derive the new, "flex-free" 
Relay, Beacon, and Monitor combinations we will start by taking specific combinations of 
the one-way data entering in each of the expressions derived in the previous section for 
the rigid-rotation case. These combinations are chosen in such a way to retain only one of 
the three noises fa, i = 1,2,3 if possible. In this way we can then implement an iterative 
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procedure based on the use of these basic combinations and of time-delay operators, to cancel 
the laser noises after dropping terms that are quadratic in L/c or linear in the accelerations. 
This iterative time-delay method, to first order in the velocity, is illustrated abstractly as 
follows. Given a function of time \I/ = ^(t), time delay by Lj is now denoted either with the 
standard comma notation or by applying the delay operator Di introduced in the previous 
section 

A* = * < = * (t - U(t)) (36) 
We then impose a second time delay Lj(t): 

DjDiV = tfjy = (t - Lj{t) - U{t - Lj(t))) 

~ V(t-Lj(t) -Li(t) + L i {t)L j ) 

~ + ifjjLiLj (37) 

A third time delay Lfc(t) gives: 

DkDjDiV = % jk = V(t - L k (t) - L,(t - L k (t)) - Li(t - L k (t) - Lj(t - L k (t)))) 

~ * >iifc + ^^fcfL^Lj + L k ) + Lj-Lfc] (38) 

and so on, recursively; each delay generates a first-order correction proportional to its rate 
of change times the sum of all delays coming after it in the subscripts. Commas have now 
been replaced with semicolons to remind us that we consider moving arrays. When the 
sum of these corrections to the terms of a data combination vanishes, the combination is 
called flex-free. 

Also, note that each delay operator, has a unique inverse, Z)r , whose expression can 
be derived by requiring that D\ = I, and neglecting quadratic and higher order velocity 
terms. Its action on a time series *&(t) is 

D^(t) ee tf(t + Li{t + Li)) . (39) 

Note that this is not like an advance operator one might expect, since it advances not by 
Li(t) but rather Li(t + Li). 

A. The Unequal- Arm Michelson 

Here we re-derive the generalized unequal-arm Michelson combination X\, by im- 
plementing our method based on the use of time-delay operators. We use again X\ as an 
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example for showing the effectiveness of this alternative and powerful method for deriving 
TDI data combinations accounting for rotation and time- dependence of the LISA arms. 

Let us consider the following two combinations of the one-way measurements entering 
into the X observable given in the previous section, evaluating them for the noises fa only 
(equation IT5j) 

+ Vl3;2) + (V21 + f?12;3');2'2] = [D 2 D 2 ,D 3 ,D 3 - I) fa , (40) 
[{r}21 + ?7l2;3') + (r}31 + »?13;2);33'] = [D 3 >D 3 D 2 D 2 > ~ I) fa . (41) 

If the time delays were constants, so the operators on the right would permute freely, simply 
differencing of equations (}4*Ul l4*Tj) eliminates fa and indeed is just X. If they do not permute, 
from equations (|4*0| l4*T|) we can use the delay technique again to write the following expression 
for X l 

X\ = \D 2 D 2 ,D 3 iD 3 - I] [(?721 + ??12;3') + (^31 + ^7l3;2);33'] 

- [D 3 'D 3 D 2 D 2 i - 1} [(r) 31 + T] 13 . 2 ) + (T] 21 + ?7l2;3');2'2] 

= [(7731 + ^13;2) + (V21 + ?7l2;3');2'2 + (V21 + ??12;3');33'2'2 + (^31 + »?13;2) ;33'33'2'2] 

- [(V21 + Vl2;3>) + (V31 + ^13;2);33' + (^31 + ??13;2);2'233' + (%1 + Vl2;3' ) ;2'22'233'] -(42) 

After substituting equations (fTUl fTTl [TJJ) into equation (j4*2^) , we obtain the final expression 
for X l Q 

X\ = [(s 3 i + Sr 3 - 2 ) + [S 2 i + Sl2;3');2'2 + ($21 + Sl2;3');33'2'2 + ( s 31 + Sl3 ;2 ) ; 33'33' 2 ' 2 ] 
— [( s 21 + Sl2;3') + {S31 + Sl3;2j;33' + («31 + Sl3;2j2'233' + ($21 + Sl2;3');2'22'233'J 

+ - [(T21 - r 3 i) - (r 2 i - r 3 i);33' - (r 2i - r 3 i) ;2 '2 + (r 2 i - r 3 i) ; 33/3 3 / 2 / 2 

+ ( r 21 — r 3l);2'22'233' ~~ ( r 21 — r 31 ) ;2'233'33'2'2] j (43) 

As usual, A 2 and A3 are obtained by cyclic permutation of the spacecraft indices. This ex- 
pression is readily shown to be laser-noise-free to first order of spacecraft separation velocities 
Lf. it is "flex-free". 

B. The Relay 

In order to derive the expressions for the generalized Relay combinations (Ux,U 2 ,U 3 ) 
valid for the realistic kinematical configuration of the LISA spacecraft, let us consider the 
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following combinations of the data that enter into the expression for U given in the previous 
section 

[l7l2 + ??31;3 + ??23;23 + ^32;123 ~ ^32] = [D 3 D 2 D 1 - I] Dy(f) 3 , (44) 
[Vl2;lV + V23;V + ^31;311'] = Dy [DiD 3 D 2 ~ I] 03 • (45) 

In each case we evaluate them for the noises 0, only, as these are what our combinations 
must remove. The expression for JJ\ is then given by the following linear combination of the 
properly delayed equations 



Ul = Dy [DiD 3 D 2 - I] [?7l2 + ??31;3 + ^23;23 + ^32;123 ~ ^32] 

- [D 3 D 2 D 1 - I] Dy [?7i2;ll' + ??23;1' + %1;311'] (46) 
= [^12;2311' + ? ?31;32311' + ^23;232311' + ^32;1232311' ~~ %2;2311'] 

— \f]Vl\V + ??31;31' + ??23;231' + ??32;1231' — ^32;l'j 

+ \J]\2;\VV + ??23;1'1' + ^31;311'1'] ~ [f?12;ll'l'123 + ^23;1'1'123 + ^31;311'1'123] (47) 

which, in terms of the one-way measurements and Ty, becomes 

Ul = [Sl2;2311' + s 31;32311' + s 23;232311' + s 32;1232311' — s 32;2311'] 

— [Sl2;l' + S31 ; 31' + S23;231' + S32;1231' ~ s 32;l'] 

+ [Sl2;ll'l' + S23;l'l' + S31;311'l'] ~ [si2;lVl'123 + S23;1'1'123 + S31;311'1'123] 

+ ^ K 7 " 32 _ r 12);2311' + ( r 21 — T3l);32311' + ( r 13 — T"23);232311' — ( r 13 — T"23);1'1232311' 

+ (Tl3 - T23);1'2311' ~ (t~32 ~ n 2 ) ; y - (t 2 i - r 31 ) ;3 i/ - (r 13 - r 23 ) ; 231' 

+ (Tl3 — T 2 3) ; l'l231' + (T32 — Tl2);ll'l' + ( T 21 ~ 7~3l);311'l' ~ (^"32 ~ TI2) ;11'1'123 

_ (Tl3 — T23);1'1'123 — ( r 21 ~ 7~3l);311'l'123] (48) 

with U2, U3 obtained by cycling the spacecraft indices. It can readily be verified using 
equations IHHJ) that the laser noise remaining in this combination vanishes to first order 
in the spacecraft relative velocities Lj. 
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C. The Beacon 

In order to derive the expression for Pi let us consider the following data combinations 
entering into the expression for P given in Section III 

[?723 + VS2(1 + Vi3-,i'i - Vi3]. 3 = D 3 [D X D V - I] £> 2 '01 , (49) 

[??32 + V23;l> + Vl2,ll> - V12} ;2 ' = D 2' [DyD 1 - 1} D 3 (j)i , (50) 

where the expressions on the right-hand-sides follow from the chosen order of the indices 
appearing on the left-hand-side of the above equation. By applying our method we obtain 
the final expression for Pi 

Pi = D 2 > [D V D 1 - I] D 3 [7/23.3 + 77 32; l3 + ?7l3;l'13 ~ Vl3;3\ 

-D 3 [DxDv - I] D 2 > [r] 32 -2' + V23;V2' + ?712,11'2' ~ Vl2;2>] (51) 
= b723;3311'2' + ^32;13311'2' + ^13;1'13311'2' - ?7l3;3311'2'] 

— [^23;332' + ??32;1332' + r 7l3;l'1332' _ ?7l3;332'] 

+ [)732;2'2'3 + ^23;1'2'2'3 + ?7l2,ll'2'2'3 ~~ ?7l2;2'2'3] 

— [^32;2'2'1'13 + ^23;1'2'2'1'13 + ^12,11'2'2'1'13 — ?7l2;2'2'l'13] (52) 

Equation (J32J) can be rewritten in terms of the one-way measurements s^j, T^j 

Pi — [S23;3311'2' + S32;13311'2' + Sl3;l'13311'2' — Si3 ; 33H'2'j 

— [S23;332' + S32;1332' + Sl3;l'1332' — Sl3;332'] 

+ [S32;2'2'3 + s 23;l'2'2'3 + S;L2,11'2'2'3 — Sl2;2'2'3] 

— [S32;2'2'1'13 + S23;1'2'2'1'13 + Sl2,ll'2'2'l'13 — Sl2;2'2'l'13] 

+ 2 [( Tl3 ~~ T 23);3311'2' — (T13 _ T23) ;1'13311'2' _ ( r 21 ~ 7~3l);2'l'13311'2' + (^21 ~ T31 ) ; 2'3311'2' 

— ( r 13 — r 23);332' + ( r 13 — r 23) ;1'1332' + ( r 21 — r 3l) ;2'1'1332' ~~ ( r 21 — r 3l);2'332' 

+ ( r 32 — r 12);ll'2'2'3 ~~ ( r 32 — r 12);2'2'3 ~~ ( r 32 — r 12);ll'2'2'l'13 + ( r 32 — r 12);2'2'l'13] (53) 

with P 2 , P3 obtained by cycling the spacecraft indices in Eq. Substituting into equation 

(153)) the laser phase noise terms entering the s^- and 7y, and applying the expansion rules 
of equations (JHol - I3HJ) . it can again be shown that, to first order in the systematic relative 
velocities of the spacecraft, laser phase noise is eliminated. 
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D. The Monitor 



The derivation of the generalized "Monitor" combinations (Ei,E2,E%) is more compli- 
cated, and rather different from the derivations shown in the previous two subsections. One 
peculiarity of these combinations is that they are not unique. It is indeed possible to derive 
different expressions for each Monitor combination. These combinations cancel the laser 
noises to the required order in the velocities, and they differ only in the number of terms - 
delayed data time series - they include. We have derived expressions with 64, 32, and 40 rj 
terms (which we do not provide here). The expression we present in this section shows the 
same number of r/-terms (16) as X x , P x , and U x . 

Let us consider the following terms entering into the expression for E derived in the previous 
section. 

%1;1'1 " Tfei = [/ - D X D V ] 01 - [/ - DiDy] Dyfa (54) 
-7732,3' - 7/23;l'3' = D 3 > [I - DyD X ] (j) 2 , (55) 
7/31 " 7731;ll' = - [/ - DyDx] <j) X + [I - D V D X ] D 2 <j) 3 , (56) 
7/23;2 + 7732;12 = ~D 2 [I ~ D X D V ] 3 (57) 

If the delay operators were constant and commuted, adding these four equations would cancel 
all laser phase noises and give E. Otherwise the above expressions can be first combined in 
pairs to remove the 2 , 3 noises in two shear- free ways 

Dy [I - DyDx) [r/21;l'l ~ T)2l] ~ [I ~ D^y] Dy [7732,3' + 7723;1'3'] 

= Dv[I-D v D x \[I-D x D^<l> x (58) 



D 2 [I - D X D V \ [?73i;ii' - 7/31] - [I - D V D X ) D 2 [r] 23 . 2 + rj 32 . X2 ] 

= D 2 [I- D X D V ] [I - D V D X ] <p x (59) 

Now we could of course repeat our iterative procedure by properly using the delay operators 
shown on the right-hand-side of equations (p)%l IBTJj) . and derive the final expression for E x . 
However, this expression would include 64 77-terms. An alternative, and more elegant way to 
derive an expression for an E x that has only 16 77-terms is by noticing that if we first apply 
inverse operators D^ 1 and D^ 1 from equation (fBl^) to both sides of equation (|5T?j) and (|5T)|) 
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respectively, and then take the difference of the resulting expressions, we get the following 
simpler expression for E\ 

E x = [I - D v Dx) [r/2i ; i'i - 5721] - Dy 1 [I - D X D V ] D 3 , [7732,3' + ^23 ; i'3'] 

+ [/ - DiDy] [7731 - 7/31.11/] + D 2 l [I - D V Dx\ D 2 [7/23;2 + r/32;l 2 ] , 
= [^31 - %L;11' - ^31;1'1 + %l;ll'l'l] ~ [V21 ~ V21;W ~ ^21;1'1 + ^21;1'111'] 
+ [V23;2 + 7 ?32;12 — V32;3' ~ ?723;1'3'] ~ [ ? ?23;2211'2 + %2;12211'2 — fe^'l'lS' 

— %3;1'3'3'1'13'] • (60) 

In equation ()60|) we have introduced a bar over some indices for representing the action 
of the corresponding inverse operator. It is easy to verify, to first order in the spacecraft 
relative velocities, that the above expression is laser and optical bench noise-free. 
Equation (|6U|) can be recast in terms of the one-way measurements s^- and 

E\ = [S31 — S3i ; n' — S3i ; i/i + S3i ; ii'i'l] — [ s 21 — S2l;ll' — S21;l'l + S21;l'lll'] 

+ [ s 23;2 + S32;12 — S 3 2;3' — S 23;1'3'] ~ [ s 23;2211'2 + s 32;12211'2 — S 32;3'3'1'13' 

— S23;1'3'3'1'13'] + ^ K 7 " 32 ~~ T ^);3' ~ ( r 32 ~ 7~12);3'1'1 ~ (^32 — 7~12);3'11' 
+ (T"32 - Tl2);3'l/lll' + (t"21 ~ r 3l) ~ (t"21 ~ T3l);ll' - (t"21 ~ r 3l);l'l 

+ ( r 21 — 73l);ll'l'l + ( r 13 — r 23);2 — ( r 13 ~~ r 23);l'12 — ( r 13 — r 23);2211'2 

+ (Tl3 - T 2 3);l'12211'2] , (61) 

with E2, E3 obtained by cycling the indices. 
E. The ( Combinations 

The expression for d derived in the previous section cancels the laser noise exactly un- 
der the assumption of constant time delays. Although perfect cancellation is no longer 
achieved when relative motion between the spacecraft is included, the ordering of the delays 
determined by our derivation of the expression for (j given in the previous section implies 
a minimization of the magnitude of the remaining laser noises at least for the equilateral 
LISA case. 

Consider the expression for (i given in Eq. (|34|). now however with semicolons rather 
then simple colons. After some algebra, it is possible to derive the leading order contribution 
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due to the residual laser noises remaining into d: 

Cl - [02,4L - 02,3L] (L 3 - Lx) L + [03,4L - </>3,3l] (Z<2 ~ L X ) L , (62) 

where we have assumed the arm lengths to differ from a nominal LISA arm length L by only 
a few percents Q ■ 

This residual laser noise can be compared with the optical path and proof mass noises 
in C\. Using the derivative theorem for Fourier transforms and taking the arm lengths to 
be the same, the spectrum of the residual laser noise in (j can be expressed in terms of the 
spectrum of the raw laser phase noise, S^, and the velocities Lf 

16 tt 2 f sin 2 (7r/L) S*(f) [{L 2 - Uf + (L 3 - L,) 2 } L 2 (63) 

From Section III and 0, the spectrum of £1 due to proof mass and optical path noises is 
equal to: 

4 sin 2 (7r/L) [24sin 2 (7r/L) S proof mass (/) + 6 S opt - path (/)] . (64) 

In Figure 2 we compare the spectrum of residual laser noise in and the optical path and 
proof mass noises in (j. The parameters used were: 30 Hz/v^Hz for the raw laser frequency 
fluctuations, 3 x 10~ 15 m/sec 2 /\/Hz for the proof mass noise, and 20 x 10~ 12 m/-\/Hz for 
aggregate optical path (shot noise, beam-pointing noise, etc.) noise. All the above spectra 
are one-sided. Figure 2 shows this comparison using nominal L = 16.67 sec. arm lengths 
and (pessimistically) taking the velocity differences to be 10 m/sec. From Figure 2, the 
residual laser noise in d for a shearing array (but with the time delays applied as given in 
equation (JMJ)) is ~ 7 dB below the optical path and proof mass noises. 
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